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SijT , Final state hyperon-nucleon interaction in strangeness photoproduction of the 

(^ I deuteron is investigated making use of the covariant reaction formalism and the 

^^ ■ P-matrix approach to the YN system. Remarkably simple analytical expression 

for the amplitude is obtained. Pronounced effects due to final state interaction are 
predicted including the manifestation of the 2.13 GeV resonance. 
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1 Introduction 

Up to now most of the information on the hyperon-nucleon (YN) inter- 
action has been obtained either from hypernuclei or from K^d and n^d 
reactions. After several decades of studies our knowledge on the YN system 
is still far from being complete. 

Recently the interest to the AA^ - SA system has flared again in con- 
nection to the expected CEBAF experimental results on the kaon pho- 
toproduction on the deuterium [1]. The final state YN interaction (FSI) 
plays an important role in the ^d -^ K^Yn, Y = A,Ti^ reaction. Therefore 
high resolution photoproduction experiments can substantially deepen our 
understanding of the YN dynamics. 

The problem of FSI in jd -^ K^Yn reaction has been addressed by 
several authors starting from the pioneering paper by Renard and Re- 
nard [2,3,4]. Two novel features differ the present work from the previous 
studies. First, covariant formalism based on direct evaluation of Feyn- 
man diagrams is used which allows to analyse the data beyond the region 
of low energy and low momentum transfer. Second, the YN interaction 
is described within the P-matrix approach which takes into account the 
subnuclear degrees of freedom and disentangle the dynamical singularities 
from kinematical threshold effects [5] . The P-matrix analysis of the YN 
interaction was presented in [6] (see also [7]) and we shall use the set of 
parameters from [6]. 

In our previous publications [8] , [9] we have announced some preliminary 
results without proving the central assertion, namely that the FSI effects 
allow remarkably simple evaluation within the P-matrix approach. In the 
present paper we give a detailed derivation of this result. Our central result 
is that YN FSI in strangeness photoproduction of the deuteron is quite 
significant and that P-matrix provides a remarkably transparent approach 
to this phenomenon. 

The organization of the paper is the following. In Section 2 we present 
kinematical relations relevant for the reaction and display equations for 
the cross sections. Analysis of the plane waves approximation is given in 
Section 3. In Section 4 we present a summary of the P-matrix approach 
to the hyperon-nucleon FSI. In Section 5 the loop diagram with YN FSI 
is calculated and expressed through the enhancement factor in the spirit 
of the Migdal- Watson FSI theory. Numerical results and calculations are 



presented in Section 6. 

2 Kinematics and Cross-section 

The reaction ^d -^ K'^Yn, y = A, S'^ is a 2 ^ 3 process. The correspond- 
ing kinematical relations as well as expressions for various cross-sections 
and distributions are well known (see e.g. [10]). However we find it apropri- 
ate to include a list of basic equations in order to facilitate comparison with 
other theoretical papers and to provide a framework for the presentation 
of the experimental results. 

We start with a general expression for the 2^3 reaction 

d- = {^-)^'\k + p,-p^-py-Pr.)^j-^^^^^^^^ (2.1) 

where k is the four-momentum of the initial photon, / is the fiux factor 
which in the deuteron rest frame is equal to kM^.., and the Dirac spinors 
are normalized according to uu = 2m. The number of essential final state 
Lorentz scalar variables for 3-body final state is 4 with the standard set 
being 

si = {pk + Prf, S2 = (py + Pn)^, ti = (k- PKf, t2 = (pd - Pnf- (2.2) 

In addition the total energy squared s = {k + Pd)'^ = {pk + Py + Pn)'^ will 
also be used in equations to follow. 

The physical quantity we are interested in is the double differential 
cross-section 

,2, ^ ^^ ^ ^^<^^^!!(f2lJ]±<l f au'yJTf . (2.3) 
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Here /c, p|^, Ek and Q^^ correspond to the deuteron rest system (LS) 
with z-axis defined by the incident photon beam direction k. The solid 
angle ily^ is defined in the Yn center-of- momentum system (CMYn). The 
quality A(x, ?/, z) is the standard kinematical function A(a:, ?/, z) = x'^ — 2{y+ 
z)x + {y — z)'^. Integration over dUyn can be replaced by integration over 
dt2dsi since linear relations hold between ^2 and cos^* and si and cos if* 
correspondingly. With such a replacement we can write d'^a in terms of 
the Chew-Low plot as 

^'^ - 1 f,„_,. I^*-"^'!.,. . (2.4) 



ds2dh 2i07r4A(s,0,Mj)V ' [-A4(ti, si, S2, ^2)]!/^' 



where A4 is a 4 x 4 symmetric Gram determinat(see [10] for an explicit 
expression) which is a quadratic in any of the variables ^1,51,52,^2- The 
region of integration in equation (2.4) has to satisfy A4 < 0. For 52 and ti 
fixed (which corresponds to fixed |pfc| and Q.k) the boundary of integration 
A4 = can be solved for 5i in terms of t2'- 

Si = ml+m^.- — {(ti+m|)(ti+m2.-t2)±Ai/2(ti,m|,0)Ai/2^t^^^2^^^2)}. 

(2.5) 
Expressed in terms of 5^ , the function A4 reads 

I6A4 = A(52,tl,M|)(5i - 5+)(5i - 5r), (2.6) 

Actual calculations presented below were performed making use of (2.4) 
since integration over invariant variables ^2 and 5i is more straightforward 
than integration over the angle ^Yn needed in (2.3). 

3 The Plane Waves Approximation 

We shall use covariant relativistic approach to calculate the amplitude T 
of the 7(i -^ K^Yn process. The amplitude will be approximated by the 
contributions of the two diagrams, namely the tree (pole, or plane waves) 
graph and the triangle graph with Yn FSI. It will be shown that within 
approximations to be specified below the sum of these diagrams reduces to 
the plane waves diagram multiplied by an enchancement factor depending 
on 52 and therefore factorizing out of the integral in (2.4) 

To calculate the plane-waves diagram two blocks have be specified: (i) 
the elementary photoproduction amplitude M^^ on the proton, and (ii) 
the deuteron vertex Td- The elementary amplitude used in the present 
calculation was derived from the tree level effective Lagrangian with the 
account of several resonances in 5, t and u chanels [11]. This amplitude has 
the following decomposition over invariant terms [12] 

M^^ = uyT. AjMj{s',t',u')up , (3.1) 

where 

s' = {k + ppf, t' = {k-pK)\ u' = {k-pYf. (3.2) 

The decomposition of the deuteron vertex function T^ in independent 
Lorentz structures has the form [13] 
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S^ipdA). (3.3) 



Here S^{pd, A) is the polarization 4- vector of the deuteron with momen- 
tum pd and polarization A. The vertex (3.3) implies that the deuteron as 
well as the spectator neuteron are on mass shell while the proton is off its 
mass shell. In the nonrelativistic limit (but with relativistic corrections) the 
functions ipi and if 2 may be expressed in terms of the deuteron S- and D- 
waves [14] 

m%. . niN. , 1 .^ TUN. . .^ .. 

^2 = wrri^y^us + ud) , (3.5) 
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with En = yp^ + mj^. The normalization condition is the following 

^f *|HI[4W + .4(p,.)] = 1. (3.6) 

Now we can write the following expression for the tree diagram 

T^'^ = UYl(x^AjMj{s',t',u')\ S{pp)rd\ul , (3.7) 

where S{pp) is the proton propagator and u'^ is a charge conjugated neutron 
spinor. The deuteron vertex in (3.7) may be substituted by the relativistic 
deuteron wave function according to [15] 

iljd=[2{2nfMd]'^'S{pp)rd. (3.8) 

Then (3.7) can be rewritten in the following schematic form 

rW = [2{27TfMd]^^^M^^iljd, (3.9) 

where il^d is the relativistic deuteron wave function discussed at length 
in [15], and where summation over magnetic quantum numbers is tacitly 
assumed. 

In the numerical calculations of the plane waves diagram use has been 
made of the Gross relativistic deuteron wave function [15]. Finally we 



mention that nonrelativistic reduction of Eqs. (2.3), (2.4) and (3.9) is 
straightforward and yields the standard impulse approximation 

^ = 6jffi!7i/'^«-|M^-nV,5t. (3.10) 



where the nonrelativistic deuteron wave function normalization condition 
is 

Jd'pnm''iPn)? = l. (3.11) 

4 Hyperon-Nucleon interaction in the P-matrix ap- 
proach 

Before we consider the triangle diagram with Yn rescattering it is worth to 
recapitulate the main results of the P-matrix approach to hyperon-nucleon 
interaction [6] . Since FSI is essential only at low relative momenta of the 
Yn pair the nonrelativistic threatment is quite legitimate. 

First we recall the interpretation of the P-matrix in terms of the under- 
lying coupled-channel quark-hadron potential [16]. Namely, the Jaffe-Low 
P-matrix description is equivalent to the coupling between hadron and 
quark channels via the nonlocal energy dependent potential of the form 

v,,, = T. ^i'^^f'\ (4.1) 

where P„ are the energies of the six-quark "primitives" [17,5] (P-matrix 
poles), and the form factors /n(^) are given by 

fn{r) = ^<5(r - 6). (4.2) 



Here h is the bag radius [17,16,5] and |c„p = Xf^/2fiYn, where A„ are the 
residues of the P-matrix [17,16,5] 

P = Po + Et7^- (4-3) 

As it was shown in [6] a single primitive at P„ = 2.34 GeV is sufficient 
for a high quality description of the existing YN experimental data, so 
that the subscript n will be omitted. To avoid lengthy expressions (4.3) is 
written for the case of a single hadron channel, namely AA^, inclusion of 



the AA^ — T,N conversion is straightforward [6] and the final results with 
the account of the UN channel will be presented below. 

The scattering problem corresponding to the interaction (4.1), (4. 2) can 
be easily solved. In particular for the momentum space half-of-shell am- 
plitude F(p,p'; E'), E' = p^/2iiYn 7^ P^/'^l^Yn and for the Jost function 
f{—p') one gets 

Fip,p;E) = --Xb-^d {E)^^. (4.4) 

d{E') = E'-En + —e'P'^ sinp'b, (4.5) 

P 

^"'(-P) = p, /'_Jx^:,, ■ ,, ■ (4-6) 

E' — En + ^e*P " sm p'b 

These results will be used in the next section devoted to the calculation of 
the triangle diagram with Yn FSI. 



5 Triangle Diagram 

Next we consider the triangle (loop) diagram with Yn FSI. The general 
expression for this graph reads 

T^'^ = / ^ Mpy) I [e AJM^ s{pp)r,cs{p^)TYnS{pY) I Mpu) ■ 

Here C is the charge conjugation matrix, Tyn is a hyperon-nucleon ver- 
tex, this vertex being "dressed" by the corresponding spinors constitutes 
the hyperon-nucleon amplitude F (see (4.4)). 

The comprehensive treatment of the triangle diagram is a purpose of 
the forthcoming publication while here we resort to some approximations 
with the aim to expose the effects of the FSI. Only positive frequency 
components are kept in the propagators S{pn) and S{py) in (5.1), while the 
propagator S{pp) together with F^; is lumped into the relativistic deuteron 
wave function i/j^ [15]. Then integration over the time component dp'^is 
performed. Thus we arrive at the following expression for T^^' 

r(') = 2M(2.)^2Afj'/7^ ^^"'^^^fa'),ffa.^-^') , (5.2) 

•^ zTT r p — p — lyj 



where p and p' are the YN cm. momenta before and after the FSI, p = 
p^ — |(k — Pk)- The quahty F{p,p'; E') is the half-off-sheU Yn amphtude 
given by (4.4)-(4.5). The arguments of the elementary photoproduction 
amphtude M"^^ in (5.2) are specified by (3.1)-(3.2) but in the kinematical 
region where FSI is essential M^^ can be considered as point-like and hence 
factored out of the integral (5.2) at the values of the arguments fixed by the 
energies and momenta of the initial and final particles (i.e. at the values 
of the arguments corresponding to the plane- wave diagram). Thus (5.2) 
takes the form 

r« = -A7i\%\2nf2Md]^'^d-\E'f^^M^''j{p'), (5.3) 






where 

dpn sin pb ipdjPr. 

(27r)3 pb p2-p'2-^o' 

In the complex p plane integral (5.4) picks up contributions from the poles 
of the propagator and from the singularities of the deuteron wave func- 
tion. The dominant contribution stems from the poles of the propagator. 
Numerical calculations using realistic deuteron wave function will be pre- 
sented elsewhere. Here we present simple arguments illustrating the above 
conclusion. Consider zero-range deuteron wave function which in the mo- 
mentum space has the form 



7r(p2+a2) 



where a = ym^lE'sl ~ 45.7 MeV. Recalling that Pn = P + Q/2, Q = 
k — px, and performing in (5.4) the angular integration, one gets 



J[P) = / TTTT^ ;— (P ~P ~ ^0) ^^In 

^^ TV Jo (27r)3 pb ^ ' pQ 



{p+^f + a' 



{p - W + a^ 



In making the approximation (5.6) use was made of the fact that FSI is 
effective at j96 < 1 while the typical momentum transfer is Q ^ 1/6 and 
Q ^ a. Formally this is reflected in the fact that the integral J{p') is 
damped at high values oip due to the formfactor sin pb/pb. the remaining 



integral (5.5) is trivial and yields 

Then using (5.3) and (4.6) one gets for the loop diagram 

^"' =^'"1-1 + 7(^(1 -^)l- (5.8) 

Here T^^^ is the plane- waves amplitude (3.9), f{—p') is the Jost function 
given by (4.6) and 

X^b sin p'b 

""^ En-E p'h 

Within few percent accuracy the term 4a^c^/Q^ can be neglected since 
Q^ ^ o? and c <C 1. The last inequality is based on the set of parameters 
fitted to the YN interaction [6], namely A^ = 0.014 GeV^ and E^ = 2.34 
GeV. For this set of parameters c ^ 0.2 in the region where FSI is impor- 
tant. 

Therefore for the sum of the plane waves and loop diagrams one gets 

7^(i) + 7^(0^rW/j(_y) (5.9) 

Thus we arrived at a remarkably simple and physically transparent result 
in the spirit of the Migdal- Watson FSI theory. This result might seem 
trivial since it can be immediately obtained by factoring out the deuteron 
wave function|^ out of the integral (5.4). However apriori such an approx- 
imation does not look reliable and has to be justified making us of the 
amplitude corresponding to the P-matrix as it was done above. Analogous 
derivation with the realistic deuteron wave function will be presented in 
the forthcoming publication. Inclusion of the AA^^ — SA conversion in the 
P-matrix approach is very simple [6]. Namely, in ( |5.!j| ) the two-channel 
Jost function has to be used. It is given by 

F-\-A. -P2) = ^, ^ ^'~xf^ x^, (5.10) 

E — En + pO l-^/ + wr^ 



where the indices 1,2 correspond to A A and EA channels and the values 
of the parameters were fitted in [6]. Then the final expression for the 

^At the value of p„ corresponding to the tree diagram with the same kinematical variables of the 
initial and final particles. 



amplitude corresponding to the sum of the plane waves and loop diagrams 
with couple-channel FSI included reads 

T^i) + T^O ^ T^t)/F{-p\, -p'^). (5.11) 

Thus using simple arguments we have splitted the description of the 
'jd -^ K^Yn reaction into two blocks: (i) the covariant calculation of the 
plane waves diagram with relativistic deuteron wave function and relativis- 
tic elementary photoproduction amplitude, and (ii) the calculation of the 
P-matrix enhancement factor with a set of parameters fitted to the existing 
data on YN interaction. The accuracy of this factorization is of the or- 
der of few percent which can be seen from model considerations presented 
above and will be confirmed by numerical calculations in the forthcoming 
publication. 

6 Results and Conclusions 

Here we present the results of the calculations performed according to 
Eqs. (2.4), (3.9), {^M and ( ^Xq) . In Fig. 1 we plot the double differen- 
tial cross-section (2.3) at |pfc| = 1.091 GeV and Qk = O'^ in a rather wide 
range of the photon energy covering the YjN threshold region with nearby 
2.13 GeV resonance. The behavior of the cross section is very instructive. 
Consider first the dased line corresponding to the plane waves approxi- 
mation. At this particular kinematics {Pk = 1.091 GeV, Q^k = 0°) the 
threshold of An production is at Ej = 1.4. GeV and according to (3.10) 
the factors that governs the near threshold cross section behavior are the 
phase space and the momentum dependence of the deuteron wave function. 
The sigma-nucleon threshold in Fig.l corresponds to E^ = 1.5 GeV. The 
T,N cusp in the cross section is practically unvisible. This fact was known 
long ago - see e.g. [18]. It is a typical situation since "In the experimental 
cross sections of nuclear reactions the cusps are rarely seen" [19]. To see 
a cusp one should accurately take into account FSI and moreover in order 
to be pronounced the cusp should be enhanced by nearby pole like 2.13 
GeV resonance close to UN threshold. Next we turn to the solid line in 
Fig. 1 with FSI included. The AA^ FSI interaction causes the cross section 
to rise more rapidly close to threshold as was noticed before in [3]. Then 
one observes a very clear signal at the SA threshold {E^ = 1.49) GeV 
caused by the pole at 2.13 GeV. We remind that according to the analysis 
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of [6,7] the 2.13 GeV resonance is not a genuine six-quark state but a P- 
matrix partner of the deuteron. Our prediction that this resonance can be 
investigated in strangeness photoproduction reaction is at odds with the 
conclusions made in [3] that FSI is important only close to An threshold 
and is completely insignificant at higher energies. The attribution of the 
pole corresponding to the 2.13 GeV resonance to a particular sheet of the 
multi-sheeted Riemann surface is a long-standing problem [6] and our hope 
is that photoproduction experiments can shed a new light on it. Concern- 
ing the structure close to UN threshold another remark has to be made. 
Namely, one may ask whether such structure can be generated simply by 
the triangle graph with constant AA^^ — T,N conversion matrix element. 
This question has been addressed in the literature for the case of a similar 
K^d -^ TT^Ap reaction [20,21], and it was shown that the singularity of 
the triangle diagram by itself can not explain the observed structure. Fi- 
nally a remark is due on another approximation tacitly made in the present 
work. This concerns hadronic corrections like pion-exchange to the basic 
P-matrix YN interaction as well as corrections due to nonorthogonality of 
the six-quark and YN ststes [16]. As was shown in [6] the present experi- 
mental data on YN interaction are not sufficient to feel these corrections. 
From the theoretical side such corrections modify the Jost function given 
by (4.6) and (|5.ii|) , in particular they smear the zero of the inverse Jost 
function at E = En. 
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Figure Caption 

Fig.l The double differential cross section as a function of the photon 
energy for px = 1.091 GeV, 6^;^ = 0°. Dashed hue corresponds to the 
plane waves approximation, solid line - with Yn FSI included. 
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